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Abstract 

A new family of nonlinear partial differential equations is presented. They represent a 
generalization of the hyperbolic Ernst equations for an Einstein- Mawxell-Weyl field in general 
relativity. A Backlund transformation for the system of equations under consideration is given, 
and their direct relation to the complete Boussinesq hierarchy of soliton equations is illustrated. 

1 Introduction 

Plane symmetric spacetimes are of fundamental importance in Einstein's theory of gravitation. 
They describe the nonlinear interaction following the collision of plane fronted gravitational waves 
solely, or coupled with other fields. In this context the basis of the vacuum Einstein equations 
is formed by the hyperbolic version of the famous Ernst equation ^ [2j . Many exact solutions of 
the Ernst equation and of its various extensions, which form the basis of the Einstein equations 
in the presence of certain fields, have been found |Sj. These solutions allow one to study the 
collisional interactions, indirectly. Assuming that the approaching waves are known, the direct 
problem is formulated by the characteristic initial value problem of the hyperbolic Ernst equations. 
For recent progress on the analysis of the preceding problem, as well as on the development of an 
effective integral representation of the solution in the interaction region, we refer to ^IHIISICI- The 
underlying structure which has lent to the Ernst equation a prominent position as a nonlinear partial 
differential equation is the large number of the internal symmetries, as was originally suggested by 
Geroch |S1 IHj . More recently in jJO] it was shown that the infinitesimal Geroch transformations 
for stationary axisymmetric spacetimes extend to a large class of equations which reside in the 
(2 + l)-dimensional Loewner system. The latter system encodes a variety of equations which from 
a modern perspective are considered to be integrable. 

One of the outstanding universal features of integrable nonlinear evolution equations is the fact 
that they are not isolated objects, but they arise always with a structured infinite dimensional 
hierarchy of mutually compatible evolution equations. The latter has been exploited in many ways, 
leading to the intimate connection between the associated hierarchy of an integrable evolution 
equation and the representation theory of affine Lie algebras (and loop groups), see e.g. [TT] 
and references therein. The point of view that one single equation in the hierarchy possesses a 
dominant character is misguided, and in fact, one should consider the infinite hierarchy of all 
equations together as being the integrable system. 

It is remarkable that there exist circumstances in which the entire hierarchy of an integrable 
evolution equation can be encoded in a single equation. A prime example of such an equation was 
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given in 12 , and it is represented by a scalar fourth order partial differential equation (PDE). 
The entire hierarchy of the Korteweg-de Vries (KdV) equations follows by systematic expansions 
from this scalar PDE. Moreover, in "J^ it was shown that a proper generalization of this PDE 
incorporates the hyperbolic Ernst equation for a Weyl neutrino field in general relativity. 

In this paper, a novel system of PDEs of fourth order in the independent variables is presented, 
which incorporates the hyperbolic Ernst equations for a source-free Maxwell field and a Weyl 
neutrino field. On the other hand this system represents in fact a consistent parameter-family of 
PDEs in the same spirit reported in J2], which encodes the entire hierarchy of Boussinesq soliton 
equations. To be more precise, the system of PDEs arises as the Euler-Lagrange equations for the 
variational problem associated with the Lagrangian density 



L = {u — v) + m{v) — + n{u) — . (1) 



-^+m{v) — +n(u) — 
The scalars quantities J-^ £ and J are given by 



~ f,UV ^ f,U1 ^ ~ f,uv ^ f,v ^ ^ ~ f,u ^ f,v > (^) 

where / = (/i,/2)- The components of /, are complex functions of {u,v) S C^, and n and 
m are complex parameter functions of the indicated arguments, unless otherwise stated. Partial 
differentiation is denoted by a comma or d, followed by the variable(s) with respect to which 
the differentiation has been performed. The Lagrangian L remains invariant, molulo L and null 
Lagrangians, under the projective group of linear fractional transformations on /, i.e. 



"31/1 + «32/2 + 033 



i = 1,2, and (aij) £ SL(3, C). If in addition the parameter functions n and m are constants, then 
L remains invariant under the following affine transformations: 

{u,v) ^ {eiu + e2, eiv + €2) , ei, £2 G C, ei / 0. (4) 

The above symmetries are called divergence symmetries of the associated variational problem, and 
every symmetry of this type is inherited as a Lie-point symmetry to the Euler-Lagrange equations, 

cf [HI- 

The integrability of the Euler-Lagrange equations associated with L, will be demonstrated by 
the existence of the usual paraphernalia of the integrable models, such as of an associated linear 
representation (Lax pair) and a Backlund transformation, in sections 2 and 4, respectively. The 
connection with the complete hierarchy of the Boussinesq soliton equations is addressed in section 
3, by exploiting the infinite symmetries of the field equations associated with the Lagrangian L. 
Section 5 deals with the reduction of the field equations to the Ernst equations of general relativity 
and the paper concludes in section 6 with discussion. 

The continuous field equations and the derivation of the associated Lax pair and Backlund 
transformation originate in a compatible parameter family, of multi-field partial difference equations 
(PAEs); the two-dimensional discrete Boussinesq system. For the relevant constructions we use an 
intrinsic property which enjoys the discrete Boussinesq system and consists of its multi- dimensional 
consistency. Recently in this property has been successfully exploited, to obtain a classification 
of one-field equations living on elementary quadrilateral faces. It can be argued that discrete 
integrable systems are equally important just as their continuous analogues, and in many respects 
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their study has led to new insights into the structures behind the more famihar continuous systems. 
The notion of Backlund transformations and their associated Bianchi permutability theorems, which 
emerged in the transformation theory of surfaces in ordinary space in the late nineteenth century, 
provide an excellent instance of such a consistency. Thus, the investigations of the interplay between 
the discrete and the continuous integrable systems acquire a certain significance, and in our view, 
the role of the second example of a generating PDE, becomes key. 



2 The field equations and their Hnear representation 

We first note that the Euler-Lagrange equations associated with the Lagrangian L, admit of a linear 
representation. This means that the field equations are a compatibility condition for the existence 
of the solution of an auxiliary, linear, overdetermined system of equations of the form 

ip^u = Mu,v;X)^ , ip^y = B{u,v;X)^ . (5) 

The complex r-vector ip and the r x r-matrices A, B are functions of n, v and a complex parameter 
A, called the spectral parameter. The compatibility condition of the overdetermined linear system 
© is 

- B^u + [A,B]=0, (6) 

which should hold identically for every value of A, where [ , ] denotes usual matrix commutator. 
Let us now consider the case where A and B are rank 1 matrices and their sole dependence on the 
spectral parameter A is given by the following rational expressions: 

_ a^ dj _ bj Cj 

u — X V — X 

i,j = 0, 1 ... r — 1. Then, equation ^ implies the existence of the parameter functions n(u), m[v)^ 
introduced by 

n{u) = di , m{v) = 5* Cj , (8) 

i = 0, 1, ... r — 1. Indices will be raised and lowered using the Kronecker delta symbol Jjj, and 
summation over repeated upper and lower indices is understood. Without loosing generality we 
may fix oq = = 1- Consequently, equation (jEJ leads to the following system of PDEs: 

(n - v) tti^-u - Cj {ai - bi) = 0, (9) 

(n - v) bi^u - V Cj {ai - bi) = 0, (10) 

{u — v) di^y + Cj di — V djCi = 0, (11) 

{u — v) Ci^u + CL^ Cj di — IP dj Ci = , (12) 

i,j = 0, 1, ... r — 1. An obvious consequence of the latter system of PDEs is the relation Ci^u = di^y, 
which implies the existence of the functions /«, such that 

Cj = fi^v di = fi^u 1 (1^) 

i = 0,1 . . . r — 1. In terms of the new variables fi and using equation (jSJ, the system of PDEs 
(jU-jni) reads 

{u - v)ai^y - [m + {aj - bj) f^^){ai - bi) = 0, (14) 

{u-v)bi^u- {n- (uj -bj) f^^){ai-bi) = 0, (15) 

{u-v)fi^uv + mfi^u-nfi^y + Cij{a^ -b>) = 0, (16) 
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where i, j = 1,2, . . . r — 1, and 

C'ij — fi,u fj,v ~l~ fi,v fj,u- (1"^) 

We note that in general, the (r — 1) x (r — l)-matrix C has rank 2. In the fohowing we consider 
the case r = 3. In this case, by ehminating the variables Oj, bi from the system H14() - ()16|) . we arrive 
at a coupled system of fourth order PDEs for the variables fi, i = 1,2. To this end, one solves 
the linear system ()16() for the differences ai — hi = Zi, in terms of the partial derivatives of the 
unknown functions /j. Substituting the result into the remaining equations (|14j) . (|15|) and using 
the compatibility condition hi^uv = (or equivalently ai^uv = CLi,vu) we get 

DuD, Zi = Du ' )+dJ L (18) 

\ u — V J \ u — V J 

i,j = 1,2, where D denotes total differentiation with respect to the indicated arguments in the 
subscripts. The fourth order system of PDEs (fTH|) . is precisely the equations obtained by varying 
the fields (/i,/2) in the Lagrangian L. 



3 Infinite symmetries and the Boussinesq hierarchy 

We next illustrate the infinite dimensional flows which are compatible with the linear representation 
©, and consequently with the field equations. This amounts to finding and representing the infinite 
symmetries of the field equations. To this end we follow the construction given in [Ifij . cf jl7j in 
connection with techniques such as the vertex operator. Having found the infinite symmetries of 
the field equations, the representation is accomplished in three steps. First, by introducing an 
infinite set of indeterminates = {x^,x'^, . . .), which are the parameters of the group associated 
with the symmetry algebra. Next, by imposing certain linear differential relations for the potentials 
tpi with respect to x^, in such a way that each of these relations is compatible with the linear pair 
©. Finally, by treating the indeterminates as additional independent variables, which leads 
naturally our considerations to the commuting compatible fiows, i.e. an infinite compatible set of 
PDEs for H{u,v,x''-). 

In order to find the symmetries of the field equations (|14|) - (|TH|) . it is necessary to introduce the 
matrix- valued function H, as follows 

i,j = 0, 1, 2. The existence of H is guaranteed by the field equations (|14|) - (|16j) . which now acquire 
the matrix form 

{v-u)H^uv + [H,u,H^^]=0. (20) 
The symmetry algebra of equation (|2U() is generated by all characteristic vector fields 

XQ = QdH + {DuQ) dn,^ + {D,Q) Oh^ + {DuD.Q) Oh^^ , (21) 

for which the matrix- valued function Q satisfies the linearized equation (|2()|) . namely 

{v - u) DuD„Q + [DuQ, H^^] + [H^u, D^Q] = . (22) 

Equation should hold for all solutions of equation (EOJ, and as such the field equations and 
their differential consequences are taken into account. In general, the symmetry characteristic Q 
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depends on the coordinates (u, f), on H and its partial derivatives up to an unspecified order, and 
on potentials of H. If Q is linear in the first partial derivatives of H, then equation H22|) delivers the 
local, or Lie-point, symmetries of equations ()2U() . These symmetries are given by the characteristic 

Qo = eii[J,H]+K) + e2iH^^ + H^,)+e3iuH^^ + vH^,) , (23) 

where J and K are complex constant matrices and Ei complex constants. It should be noticed 
that if n and m are parameter functions of u and v respectively, then only the symmetry with 
£1 7^ and £2 = ^3 = 0, is compatible with the symmetries of the field equations H14 p - H16p . since 
on the contrary the characteristic corresponds to the affine base transformations ^ . The following 
observation allows one to find in principle, an infinite number of symmetries of equations ()2U() . 
starting from the characteristic Qq with £1 = 1,£2 = £3 = 0. 
Equation 1)221) . can be written in the form 

{uDuQ + [Q, = Du {vD^Q + [Q, H^^]) , (24) 

which implies the existence of the matrix-valued potential Q introduced by 

DuQ = uDuQ + [Q, H^u] , A.Q = vD^Q + [Q, i/,^]. (25) 

By direct calculations and using the Jacobi identity, one finds that the potential Q is also a sym- 
metry of equations ^E^ . whenever Q and H satisfy and (|^ . respectively. In other words, 
equations H25|) define a recursion operator IZ, i.e. an operator with the property that whenever Q 
is a characteristic of a symmetry of equations ()2n)) . so is Q = IZQ. It follows from equations ()25|) . 
that we may formally write TZ as the operator 

7^ = • (ctD - D + add//) . (26) 

Here, D and d denote total and ordinary differentials respectively, * is the two dimensional Hodge 
duality operator which acts on the basis of one forms by 

*du = du, *dv = —dv , (27) 

ad^(i?) = [i?,y4], and the functions q and a are given by 

a = \{v + u), Q=\{v-u). (28) 

From the above discussion it is clear that starting with specific matrices J and K and using 
the recursion relations 

Qp+l=7^Qp, Qq = [J,H\+K, pGN, (29) 

we may generate an infinite number of symmetries of equations ()2U|) . Let = {xP , x^, . . .), p G N, 
denote the parameters of the infinitesimal symmetry transformations which are generated by the 
vector fields Xq^, i.e. 

d^pH = Qp, peN. (30) 

Since for every fixed p, the characteristic Qp defines a symmetry of equation ()2U() . the pair of PDEs 
which is formed by taking each of equations ()3()j) together with equation ()2())) . is a compatible 
system of PDEs for H{u,v, x'^) (p fixed). By virtue of this fact, every such pair of PDEs, admits 
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also a linear representation. The latter is given by the linear equations (jSj) for ip, and which for 
convenience we rewrite below 



{u - X)i/j^u = H^u'fp , {v - X)ip^^ = H^^ip , (31) 
together with each of the following linear equations 

tx^ = J^, (32) 

= -Xlp^xP - Qptp , peN. (33) 

By induction, we may easily prove that the compatibility conditions dxvduip = dudxPip and dxpdyip = 
dvdxPip are satisfied for every p £ N, whenever H satisfies equation (|2fl|) and Qp are the symmetry 
characteristics of the previous construction. 

Similarly, starting with some other matrices J and K, one obtains a different infinite hierarchy of 
symmeries of equation (|20|) . Let y'^ = (y^,y^, . . q £ N denote the parameters of the corresponding 
symmetry transformations generated by the vector fields . In analogy with the preceding 

Wq 

■ i ; - 1 "I : : ■ j': : I . thc linear equations which can be imposed in addition to the linear equations H31|) 
are 

V = Jip, (34) 

Jp^yq+l = -X'lp,y1 - Qgtp , gCN, (35) 

We now assume that the potentials H and ip depend in addition to {u,v), on and also on y'^, 
Vp, q £ N. Since and y'^ will be treated as additional independent variables, it is necessary that 

dxpdyqF - dyqdxpF = , V^, Q G N . (36) 

holds, for all functions F{u, v, x^, y'^). The aim now is to generate the compatible set of PDEs for H, 
in terms of the higher time variables x^, y'^. This infinite set is obtained by taking the compatibility 
conditions dxpdyqtp = dyqdxPip, on the linear representations (1^5]) and ((311), (ESJ- In order to 
give an explicit description and for purposes of calculation we proceed as follows. 

The leading equation (|3H|) is valid for every solution H of equation (|2()|) , whenever the starting 
matrices satisfy 

[J, J] = , and [J, K] = [J, K] . (37) 
By virtue of relations ()37() . we fix the values of the starting matrices as follows: 

) -1 ( 

{J^j) = I I , (K^j) = I -1 I , (38) 

) ( 

) -: 

{J^J) = I 1 I , (K^j) = I I , (39) 

) ( 

i,j = 0,1,2. Next, we factor out the linear set of PDEs for H, which is obtained from the 
compatibility conditions dxpdyqip = dyqdxPip with p = 0,q > 1, {q = 0,p > 1). Finally, we combine 
the remaining linear representations (|33|) and (|35j) into one, by introducing (i) auxiliary higher time 
variables 

= {w^ =x\w'^ =y\w^,w^,w'^ ...), ieZ\3Z, i>l, (40) 
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and (ii) the matrix 

A=-XJ-K. (41) 
From this starting point, and noticing that 

A'^ = -XJ-K, (42) 

the hnear infinite set of differential relations for ip, which can be simultaneously imposed on the 
linear representation is given by 

V',^! = (A +[A',H])ij, (43) 
= (A2+ [iA^)',H])^p, (44) 
V'.io'+a = — AV'^^i —H^^itp, i G Z \ 3Z, i>l, (45) 

where prime in equations H43() . (|44() denotes differentiation with respect to A. The above differential 
relations generate an infinite set of PDEs for H, with the Boussinesq equation being the first 
member. Specifically, the compatibility condition d^2d^iTp = d^idy^2ip on equations 
yields a set of six, first order PDEs for the components of H. After a simple differentiation and 
elimination process applied to the latter system, and using (fT^ to identify the main field variables 
(/i, 72), we arrive at the following system of PDEs: 

fl,xy = 2/2,yy + {f2f2,x),y — ^f2,xf2,xx — \ f2,xxy + jf2,xxxx , (46) 
fl,xx = {f2f2,x),x + \{f2,xy — f2,xxx) , (47) 

where x = w^,y = . Finally, taking the compatibility condition fi^xyx = fi,xxy on equations 
(|?7|) we obtain 

3B,yy + B^xxxx + Q{B^),xx = , (48) 

where B = —f2,x- Equation ()48|) is known as the Boussinesq equation. Similarly, higher members 
of the Boussinesq hierarchy can be derived from the system (|43|) -(|45 |) . in a systematic way. 



4 A Backlund auto-transformation 

The discussion in the present section relates to matters that had their origin in the investigations 
by Backlund near the end of the nineteenth century, concerning the transformation properties of 
simultaneous partial differential equations and the transformation theory of surfaces in ordinary 
space. There are many methods developed over the past century to obtain the Backlund transfor- 
mations for certain classes of partial differential equations, see jl8l 1191 12U1 121| for recent advances 
on these topics. However, for the construction of a Backlund auto-transformation for the relevant 
PDEs of our discussion, i.e. a transformation within the space of solutions of PDEs (|14() - (|16|) . we 
find the connection with the lattice Boussinesq equation, as a consequence of its three dimensional 
consistency, as the most powerfull. 

An upshot of recent investigation of the equations H18|) , and their compatible discrete counter- 
parts, has been the insight that both systems arise as parameter families of mutually consistent 
equations. In turn, these equations form multi-variable and multi-parameter systems, in which 
parameters and variables play a dual (even interchangeable) role, and the whole infinite system 
is characterized by an intrinsic consistency. In the light of the close relation of equation H18|) 
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with the Ernst equations, as will be explained in the next section, and since solution generating 
techniques such as Backlund transformations, have proven to be invaluable methods in finding 
physically significant exact solutions of the Ernst equations |22| I23| I24j . we now explain how a 
Backlund auto-transformation for equations H18|) can be constructed straightforward, by exploiting 
the consistency approach. 

The construction is based on the compatible discrete system as this was derived in . In order 
to make contact with the latter system, it is necessary to retract for the moment the dependence of 
the parameter functions n, m on u, v, respectively. ^ Let us now consider an elementary plaquette 
of a two-dimensional lattice, as shown in figure 1. The dependent variables (fields) are assigned 
on the vertices at sites (n,m), and (u, f) which now play the role of the lattice parameters, are 
assigned on the edges. The updates of a lattice variable F along a shift in the n and m direction 
of the lattice will be denoted by Fjij , F[2] respectively, i.e. 

=F(n + l,m), Fpj = F(n, m + 1), F[i2] = F(n + 1, m + 1) . (49) 

The partial difference equations which are compatible with the partial differential equations H14|l - 

F , 



F 



[2] 



-© F 



[12] 



Figure 1: An elementary quadrilateral for the discrete Boussinesq system. 

H16|) . in the sense that these two systems share a common set of solutions, is given by the following 
parameter family of algebraic equations: 

(50) 
(51) 

(52) 



^[1] = 




V, 




UU^2]- 


V, 


w = 




- V[12] + 



u — V 



[2] 



The system of equations (|50p - (|52p . relates the values of the fields assigned on the four vertices of 
the quadrilateral, where now F = (U, V, W), and is called the discrete Boussinesq system. For the 
Cauchy problem on a staircase, we should impose on the black points of the elementary plaquette 
of figure 1, the initial values {U,V,W), (C^[i], Vjij), (C^p], Vpj) only, since from these data and the 
discrete equations (f5n|) - (|52|) . the values (VF[i], W[2]) and (?7[i2], ^12]) W^[i2]) are determined uniquely. 
In particular, using equations (fSO]) . (|^T|) we find that 



V, 



12 



[1] 



[2] 



2] 



(53) 



^Relative to the continuous field equations, this restriction is not essential since the Backlund transformation that 
we derive below, is also valid with n, m parameter functions. 
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and subsequently the values VF[i2] and V[i2] can be found uniquely, from equation (ISU]) (or (|^T|) ') 
and equation (|52|) . respectively. The fields which are assigned on the vertices with black points, 
are identified with the fields of the continuous PDEs, as follows: 

W = fi, U =-f2, V = -Hu , (54) 
C/[i] = ai, U[2] = h V[i] = 02, V[2] = 62 • (55) 

In j25j it is shown that the discrete Boussinesq system represents a consistent parameter family 
of two-dimensional lattice equations, in the sense that equations (|5()|) - (|52|) can be imposed on all 
faces of a three-dimensional lattice in a consistent way. The third direction of the generated lattice 
has a dual role. On the discrete level, is taken as auxiliary (spectral) and the system is linearized 
in the corresponding variables i7[3], Vjs] and their shifts. On the other hand, the discrete equations 
which relate the fields and their shifts in the third direction, define a map in the solution space of 
the continuous field equations, by virtue of the identifications (|51)) - (|55() . 

Specifically, extending equations (|50)) - (|52|) into a third direction associated with a new lattice 
parameter (identified here with the spectral parameter A) we get 

Wis] = UU[3]-V, (56) 

f/[13] = I , ^ = t^ ^[13] -^13] + 77^^^, (57) 





-V 


^[11- 






Urn 




^3] 




Urn 



V[i3] 


+ 


U 


- A 






V[23] 


+ 


V 


- A 


Urn 





^[2] - t^[3] ^[3] - ^[2] 

The auxiliary variables C/[3],Vi3] are identified as projective variables and their relation with the 
spectral variables tpo,tpi,tp2 is given by 

By straightforward calculations one easily verifies that the linear system © for the potentials 
("001 ■01) V'2)j implies the Riccati system for the potentials {qi,q2) = (f^[3]) ^[3])) 

(u - X)qi^u = {ai - Qi) {n - {aj - qj)f^u) > (60) 
{v - X)qi^^ = {hi - qi) (m - {hj - qj)f^^) , (61) 

i = 1,2. In virtue of the identifications H54() - H55() . we now interpret the shift of the fields in the 
third direction, as a map ip in the solution space of equations (fTl|) - (|16j) . Explicitly, using equations 
((nZI)-(EHl) and the identifications (|51)) - H55() . the map ip is defined by 

M)=Hi2-f2qi, V{h) = -qi, (62) 
ip{ai) = — — — , ip{a2) = -fi - f2 v{ai) H , (63) 

oi -qi qi- ai 

V{h) = , ^{h2) = -h - f2 ^{bi) + . (64) 

hi -qi qi- hi 

Restricting now our considerations on the continuum level, the above equations define a Backlund 
auto-transformation, i.e. a transformation that generates new solutions of the field equations from 
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other known ones. Below we summarize the resulting solution generating procedure, which can be 
verified by lenghy but straightforward calculations. 

Let (/i,/2) be a known solution of equations (fT8|) . From this solution and using equations 
((TH) - (fT^ . one determines by quadrature the corresponding vector components (ai, 02) and {bi, 62)) 
up to constants of integration. Successively, the auxiliary complex potentials H12 and (^1,^2) 
can be found by quadrature from equation (fTU)) . and by integrating the Riccati system (|^ - (|61() . 
respectively. Then, a new solution (0,1,0,2), (&i, ^2) and (/i, /2) of the system (|ll|) - (|16|) . and a fortiori 
for the equations (|TH|) . is given by the algebraic relations (jf)^ - (|64j) . where fi = ^p{fi), Oi = ^{oi), 
hi = (pih), 1 = 1,2. 



5 Reduction to the Ernst equations 

We next show that the main second order subsystem of PDEs which is incorporated into the 
fourth order system of equations (|18|) . is the hyperbolic Ernst equations for an Einstein-Maxwell- 
Weyl field. In order to make this connection more transparent, we introduce auxiliary complex 
potentials {hi, /12) by the relations 

, 02-62 , 6102-0162 

hi = 7- , ri2 = r . (65) 

oi — bi oi — 61 

and which correspond to a change of the variables {02, 62), of the form 

02 = hiOi + h2, 62 = hi 61 + /i2 . (66) 
It follows from the system of equations (|14() . H15|) that 

a2,v = hiai^y, , b2,u = hibi^u, (67) 
which in virtue of equations (|66j) . yield 

h2,v = -hi^v ai , ^2,n = -hi^u 61 . (68) 

From this starting point, and eliminating the variable fi from the field equations ()14() -H16 |) . we may 
recast the latter system as a coulped set of four PDEs of second order in terms of the variables 
(oi, bi,hi, f2). This system of PDEs is given by 

Ol^uai^v m n+1 hiuf2,v , , .2 /ar,\ 

ai,uv = 2 oi^u ai.i-H (ai - 61) , (69) 

ai — 61 u — V u — V ' u — V 

, n bi^ubi,v m + l n hi^J^u, , ^2 /vn\ 

bi,uv = 2-^ H 6i,„ H 61,^, ■ ^(oi - 61) , (70) 

61 — ai u — V u — V u — V 

f2,uv = J— (6i.m/2,d — ai.vf2,u) , (71) 

ai - 61 ' 

hi,uv = 7-{bi.vhi.u — ai.uhi.v) , (72) 

ai - 61 ' ' ' ' 

where the first three PDEs follow from the system (|14() - H16|) . and equation ()72|) follows from the 
compatibility condition h2,uv = h2,vu on equations (|BH|) . 

Let now u, v be real valued, and * denote complex conjugation. It follows from the system of 
PDEs (jM |) -(f7 ^ . that the latter admits a compatible reduction by applying the reality conditions 

ai = -b\ , hi = , (73) 
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whenever the parameter functions n, m satisfy 



n + n* + 1 = 0, m + m* + l = 0. (74) 

The resulting set of equations H69() and (jTlj) . which is formed by applying the above conditions, 
are the Ernst equations for an Einstein-Maxwell- Weyl field. A more standard form of the latter 
equations is given by a dual system, which is obtained here by rewriting the system ()14() -()16 p in 
terms of the above reality conditions. These conditions imply the following constraints on the 
variables (oj, 6^): 

, V — u [v — u)f2 

Jl + Jl + /2/2 Jl + /l + J2/2 

By virtue of (|75j) . one readily notices that the system of equations (|14j) - (|16|l . decouples for the 
variables (/i,/2), leading to a second order system. In order to make contact with the notation 
employed in general relativity, we relabel the main complex variables as (/i,/2) = (.E, $). Using 
equation ((7S|) and adopting a coordinate free setting, the relevant second order system of PDEs is 
given by equations (fTHj) . and reads 



{E + E* + m''){d{Q*<lE)-if3hdE) = 2^(d£; + $M$) A *dS, (76) 
(^ + E* + $$*)(d(^*d$)-i/3Ad$) = 2^(d^ + $M$) A *d$ . (77) 

Here, i = and the one-form (3 is given by 

(3 = v{u)du + iJL{v)(iv , (78) 

where the real parameter functions the /i, are introduced by recasting the parameter constraints 
(|74|1 in the form 

n = —\ + \v{u) ^ m = — i + i/i(u). (79) 

In the context of general relativity, equations H76|) , ()77() are known as the Ernst equations for an 
Einstein-Maxwell- Weyl field [^Hl- It should be clear from the preceding discussion that, all solutions 
of the Ernst equations are embedded into the solution space of the Euler-Lagrange equations (fTH|) , 
while the converse does not hold in general. Therefore, equations (|18jl . or equivalently the auxiliary 
system ((TH) -(P3 |) . represent a generalization of the Ernst equations (f7S)) . ((77|) . 

At this point, it is worth noticing that after applying the conditions (|73p . (|79|) and their conse- 
quences to the linear representation ©, the latter reduces to the one introduced in l?71 05] . 

6 Discussion 

Apart from the Ernst equations, which possess a prominent position among the subsystems included 
into the solution space of equations (|18|) . there are also various interesting subsystems which are 
tied in with the latter set of equations. Thereof, we point out that equations H18|) . incorporate 
also the PDE which encodes the KdV hierarchy. This is accomplished by considering the case, 
= = /, which leads to a scalar fourth order PDE for the variable /. The latter equation is 
the Euler-Lagrange equations for the Lagrangian density 

r I \if,uv)'^ , 1 / 2f,u , 2 A fon\ 

L= [u- v)—^ \ m + n -T^ . (80) 

I,uJ,v u-v \ 
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In jl3j it was shown that in the case where n and m are constants, the corresponding PDE for 
/, admits similarity solutions built from solutions of the full Painleve VI (PVI) equation, i.e. with 
all parameters free. Alternatively, full PVI can be shown to arise naturally as symmetry reduction 
from the anti-self-dual Yang- Mills equations [22] , or from the Ernst equation for an Einstein- Weyl 
field 30, . Taking into account the richness of equations H18() , it is even more interesting to consider 
their similarity reductions. Such a reduction leads to a six-parameter, second order, coupled system 
of ordinary differential equations 25 , which we expect to be associated with the hierarchies of the 
PVI equation |3T] . 
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